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ABSTRACT 

We develop the improved ladder approximation to QCD in order to apply it 
to the heavy quark mesons. The resulting Bethe-Salpeter equation is expanded in 
powers of the inverse heavy quark mass 1/M, and is shown to be consistent with 
the heavy quark spin symmetry. We calculate numerically the universal leading 
order BS amplitude for heavy pseudoscalar and vector mesons, and use this to 
evaluate the Isgur-Wise function and the decay constant Fg. The resulting Isgur- 
Wise function predicts a large charge radius, p 2 = 1.8 — 2.0, which when fitted to 
the ARGUS data corresponds to the value V c \, = .044 — .050 for the Kobayashi- 
Maskawa matrix element. 



1. Introduction 



Recently there has been a great deal of interest in using B meson physics 
to test the detailed structure of the Standard Model and to reveal new physics. 
When we consider systems containing a heavy quark (such as the b quark and, 
possibly, the c quark) a new symmetry, the heavy quark spin- flavour symmetry' 11 
appears. This symmetry, which becomes exact in the limit when the heavy quark 
mass M goes to infinity, may be used to predict various properties of the heavy 
mesons; these have been systematically studied using the Heavy Quark Effective 
Theory' 2 ' 3 ' 41 For instance, when M — > oo, every semi-leptonic form factor can be 

[2] 

expressed in terms of a single universal function, the Isgur-Wise function. At the 
moment there is a great deal of work being done to understand the semi-leptonic 
weak decay process B — > D^*Hu and thence to extract the Kobayashi-Maskawa 
matrix element V c i,: [o] it is therefore important to have a reliable functional form 
for the Isgur-Wise function. Unfortunately the heavy quark symmetry tells us 
nothing about the Isgur-Wise function away from the kinematical end point. In 
order to calculate it at any other point we have to know the details of the strong 
interaction physics of the light degrees of freedom in the heavy mesons. 

In this paper we develop the improved ladder approximation to the QCD Bethe- 
Salpeter equation as an expansion in powers of 1/M. We show that it is indeed 
consistent with the heavy quark symmetry, and apply the results to calculate the 
Isgur-Wise function. In previous papers' 6 ' 71 we have shown that the improved ladder 
approximation is capable of giving rather a good description of the physics of light 
mesons, and we therefore expect it to provide a reasonable approximation to the 
physics of the light degrees of freedom in heavy quark mesons. 

We calculate the Isgur-Wise function as a one loop diagram containing two 
meson BS amplitudes and a single current insertion. The main problem which 
we must handle is the following: whilst the overall boundstate momenta must 
remain timelike, we have to Wick-rotate the loop momentum to evaluate the in- 
tegral. Then the combinations of the loop momentum and the overall boundstate 
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momentum which appear in the arguments of the quark mass functions (and also, 
depending on the improvement scheme, in the argument of the running coupling 
constant) becomes complex (non-real). Moreover, when we calculate the Isgur- 
Wise function, we work in the rest frame of one of the two mesons but the other 
must be boosted to a finite velocity, and this corresponds to the BS amplitude with 
complex relative momentum. We are thus forced to calculate the BS amplitude 
with complex arguments. We propose a new method to do this analytic contin- 
uation using the BS equation itself. Because of this complication, the numerical 
work becomes quite hard — the Isgur-Wise function, for instance, appears as a 
five-dimensional integral. 

In order to illustrate the main points of the method and obtain numerical 
results without being too much involved in such complications, we here ignore the 
running of the quark masses and work with a fixed value of the light quark mass m. 
We determine m by solving the equation m = E(am 2 ) where £ is the light-quark 
mass function determined from the ladder Schwinger- Dyson equation, and a is an 
unknown parameter. We have calculated the Isgur-Wise function for the two cases 
which we regard as typical, a = 1 and a = 4. A full treatment including the quark 
mass functions derived from the SD equation will be reserved for a later work. We 
actually find that the general shape of the Isgur-Wise function and the value of 
the charge radius p 2 are rather independent of the value of the fixed light-quark 
mass. Consequently we expect that the 'true' values obtained by including the 
mass function will not differ significantly from our present results. The heavy 
meson decay constants, on the other hand, which we also calculate do not show 
this independence of the value of m, and so we have to wait for the full treatment. 

The rest of this paper is organised as follows. In section 2 we present our 
formalism for treating the heavy meson BS equation as an expansion in powers of 
inverse heavy quark mass 1/M, and we derive the leading order BS equation in a 
component form. An exact expression for the decay constant Fb in terms of the BS 
amplitude is given in section 3. In section 4 we explain how to calculate the Isgur- 
Wise function and show that our approximation is actually consistent with current 
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conservation in the equal heavy quark mass case. In section 5 we demonstrate 
that the leading order BS equation satisfies the heavy quark spin symmetry and 
discuss the implication for the form of BS amplitudes for pseudoscalar, vector, 
scalar and axial-vector mesons. Section 6 is devoted to the details of the numerical 
calculations. In section 6.1 we fix a suitable analytic functional form for the running 
coupling constant. In section 6.2 we use this coupling and solve the SD equation 
for the light quark mass function E(x). We need S both in order to fix the light 
quark mass which we use in the BS equation, and also to fix the overall mass scale 
in MeV. We fix the scale by calculating the Pagels-Stokar pion decay constant 
and defining the result to 93V2MeV. In section 6.3 we solve the BS equation and 
calculate the decay constant Fb and the energy eigenvalues of the ground state 
and the first excited state. Finally, in section 6.4, we present the results of our 
calculation of the Isgur-Wise function £(t). 

2. BS Amplitude and Notation 

The BS amplitude of a meson boundstate \B(q)) of total momentum q con- 
taining a heavy quark \I/ and a light anti-quark ip is defined by 

<0| T* ia (x)^(y) |S(q)> = e^ X J ^ 5f Xa p( P ; q) , (2.1) 

where i,j and a,f3 denote colour and Dirac spinor indices respectively. The coor- 
dinate and the relative coordinate r M are defined by 

X» = (x^ + riyV , r^ = x^-y^ . 

where ( and rj are real parameters satisfying ( + rj = 1. p is the therefore the 
relative momentum. Although in the relativistic case any value may in principle 
be chosen for (, we make the following 'natural' choice 

M m 
M + m' 11 ~ M + m ' 

with M and m being the masses of the heavy and light quarks, respectively. With 
this choice, is the centre-of-mass coordinate of the system; it has the technical 
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advantage of making it legitimate to perform a Wick rotation in the BS equation. 
For later use we here define the 4-velocity v M of the boundstate: 

v v = /Mb , v 2 = 1 , (2.2) 

where Mb is the boundstate mass, q 2 = M B . 

It is convenient to introduce a ket notation to denote the bispinor BS amplitude: 

«X» — 4x«p(p;q). (2.3) 
The conjugate BS amplitude which is denoted by a bra state is defined by 

«X« <— > Sixap(p;q) = Si^[ Xa p(p*;q)H. (2.4) 

Note that the complex conjugate p* of the argument p appears in the RHS: this 
is necessary for the general case where the relative momentum p becomes complex 
after Wick rotation. In fact, even when p is real, the ie prescription in the Feynman 
propagator implies that the Minkowskian theory should be understood as that 
obtained by analytic continuation of the Euclidian theory. Considering the limit 
Argp — > 0+, one can see that x defined this way does indeed coincide with the 
conventional conjugate BS amplitude:* 

(S(q)| T^(y)¥ a (x) |0> = e* x J ^ & 5) X( 3 a ( P ; q) . (2.5) 
We also introduce a (non-positive-definite) inner product: 

mx)) = N c J T^tr[i>(p- q )x(p;q)} , (2.6) 

where N c = 3 comes from the trace over the colour index and the trace in the RHS 
is taken over the bispinor indices alone. 

* One can also confirm this fact by seeing that the BS equation for the conventional conjugate 
BS amplitude is indeed satisfied by our conjugate amplitude (2.4). 
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The bispinor BS amplitude \ may be expanded into invariant component am- 
plitudes. Here we are mainly interested in the pseudoscalar boundstate for which 
here there are four such amplitudes: in the present context of the heavy quark 
system, it is convenient to use the form 

X (p;q) = A+(A(u,x) + B(u,x)j)j 5 + A-(C(u,x)+D(u,x)f)'y 5 , (2.7) 

where the projection operator A± and the 4-momentum p^, the spatial projection 
of are defined by 

A± = ^ , p" = p" - (p • «y , 

and the variables u and x are 

iu — p ■ v , x = \J—p 2 = \J — p 2 — u 2 . (2.8) 

The components of the conjugate BS amplitude x are defined similarly: 

Xip) ^ -!5{A(u,x) + B(u,x)p)A+ - l5 (C(u,x) + D(u,x)p)A- . (2.9) 

Note that iu is the time component p° of p^ in the restframe of the boundstate 
where t> M = (1, 0). In this frame, u becomes real after Wick rotation, p^ — > p^, = 
(u, p), and then eq.(2.4) gives 

Xip) = -lh{A*(-u,x) + B*{-u,x)0)A+ - l5 (C*{-u,x) +D*{-u, x)^)A_ . 

(2.10) 

In view of the BS equation (given shortly in (2.13)), we can easily see the relations 
X*(-u,x) 

= X(u, x) for X = A,B, C, D, provided that the boundstate mass Mb, determined 
as a eigenvalue of the BS equation, is real. This is because when = Mbv^ is 
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real, in the Wick rotated form of the BS equation (in which the integration fac- 
tor J d 4 k/i(27i) 4: becomes real), the imaginary quantity appearing is only through 
p° = iu (or ko = iv) so that the complex conjugation of the whole equation is 
equivalent to changing the sign of u (and v). So, for the true boundstate solutions 
corresponding to real mass eigenvalue M#, we have simply 

A(u,x) — A(u,x) B(u,x) = B{u,x) 

(2.11) 

C(u, x) = C(u,x) D(u,x) = D(u,x) . 

To simplify the discussion, we use the following abbreviation for momentum 
integration and for gamma matrix multiplication on the bispinor x' 

d A p f o? 4 pe f x 2 dx du dcos 9 



i{2nY J (2tt) 4 J 8vr3 ( 2 .12) 



In this notation the BS equation becomes 

Sn( P + (q) ®S L (p- m ) ||x)) = K\\ x )) (2.13) 

where Sh and S*l are {i times) the inverse propagators of the heavy quark and the 
light quark respectively 

S R = i> + Ci - M , S L ee $ - r]4 - m , (2.14) 

and K stands for the Bethe-Salpeter kernel. In the improved ladder approximation 
K is given by the following one-gluon-exchange form: 



K\\ x ))= [ g 2 C 2 D^(p-k)^ X (k)^ , 
Jk 

g^v ~ k^k v j k 



(2.15) 



D^(k) 



k 2 



where C<i is the second Casimir given by C2 = (N 2 — 1)/2N C and g 2 is a suitable 
running coupling constant. Since the running of the coupling is mainly caused by 



7 



the gluon self-energy corrections, probably the best choice of the argument of the 
running coupling constant is to use the gluon momentum and take g 2 ((pE — ^e) 2 )> 
where p^ is Euclidean momentum; p| = — p 2 . In this paper we take for simplicity 
9 2 {Pe + ^e)- ^ e have discussed this choice of argument before 191 and have shown 
that it does not differ from the ideal case very much. The precise functional form 
of our running coupling constant will be given later. 

In the context of the heavy quark system, it is convenient to adopt the following 
boundstate normalisation 

<i?(q)|5(q / )> = (27r) 3 2t;o ( 5 3 (q-q / ). (2.16) 

This differs from the usual invariant normalisation in that we have used v in 
place of g° = Mbv°: the usual invariantly normalised state is therefore given by 
|-B(q)) \[Mb- By considering the contribution of the intermediate boundstate to 
the two heavy two light quark Greens function in the standard way, we can show 
that our BS amplitude may be normalised by requiring 

C«xiiTM ® Sl«X» - v((x\\Sn ® 7mIIx» = . (2.17) 



The BS equation gives an eigenvalue equation for the boundstate mass Mb or, 
equivalently, for the binding energy E defined as 

M B = M + m-E = CHM-CE) . (2.18) 



We now expand all quantities as power series in m/M. The light quark mass 
m (which is the constituent quark mass and comes almost entirely from dynamical 
chiral symmetry breaking) is of order Aqce>, and we therefore regard Aqce> as 0(1) 
in this expansion. Thanks to our choice of ( and i], the relative momenta p and 



8 



k are also of the order of Aqcd, and therefore also 0(1). The binding energy and 
BS amplitude are expanded as 

C £ = ^ + (_)^ + (_)^ + ..., 
x = xo + (^)xi+ (^) X2 + --- , 

with corresponding invariant amplitudes: 

Xo(p) = A+(A (u,x) + Bo(u,x)tf)7 5 + A_ (C (u,x) + D (u,x)^)-f 5 , 

(2.20) 

Xiip) = A+(A 1 (u,x) + Bi(u,x)tf)i 5 + A_(Ci(u,x) + Di(u,a;)ji)75 , 

and so on. The expansion of the quark inverse propagators is (notice that Sh starts 
with a —1-st order term) 

Sh = (-^r) ^H-l + SrO + ij^) Sm ^ ' 

Su-i = m(i - 1) = -2mA_ , 
5ho = ^ - E^i , 

m2 ( 2 . 21 ) 

Sl = S L0 + {jj)S L l + (^) 5L2 + • • • , 

S w = $ - m(l + y) = ^ - 2mA + , 

511 = ^ , 

512 = Eii/i . 

Substituting these expansions into the BS equation (2.13), we find that at —1-st 
order 

Sh-i®Slo|xo)) = , (2.22) 

at zeroth order 

Sn-i ® (^loIIxi)) + 5 L i||xo») + ^ho ® 5 LO |xo)) = ^||xo» , (2.23) 

and so on. Since Sh-i oc A_, the —1-st order equation (2.22) simply implies 
that xo oc A + so that we find Cq(u,x) = Dq(u,x) = 0. This then causes the 
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term Sh-i <8> Sli||xo)) i n the zeroth order equation to vanish. Moreover in the 
zeroth order equation, the first term, Sh-i ® <Slo||xi))j contributes only to the 
piece proportional to A_, so that we have 

A+[S m ®S w -K}\\xo} = . (2.24) 

Projecting out the A+ pieces by taking itr[7 5 A + x (2.24)] and \ tr^A+x (2.24)], 
the zeroth (or leading) order BS equation (2.24) becomes 

where the arguments p and k stand for (u, x) and (v, y) respectively, M. is a 'metric' 
matrix given by 

(iu —x 2 
—x x (iu — 2m) 

and the kernel Ko(p,k) is obtained after integrating K over the angular variable 
cos# = k ■ p/|k||p|. The explicit matrix expression for Ko(p,k) is given in the 
Appendix. 

In the leading order of the m/M expansion the normalisation condition (2.17) 
becomes 

\(xoW ® S LO \\ X o)) = 1 , (2.27) 
which when rewritten in terms of the invariant amplitudes gives 

with the metric matrix M. defined in eq.(2.26). 
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The BS equation (2.25) is a linear eigenvalue equation A ||x)) = EqM. with 
A = g 2 C2Ko + iuJ\A and metric A4. Further, both A and the metric M. are hermi- 
tian with respect to the inner product ({?p\\x} defined in (2.6): [Aij(—u, x; —v, y)}* =| 
Aj,i{y , y; u, x) = 1,2) and so on. As suggested by the normalisation condition 
(2.28), the natural inner product to this BS system is given by ((Vi-Mflx)) an d the 
norm by ||x|| 2 = ((x||-M|x)) accordingly. Indeed, the BS equation (2.25), with the 
help of the hermiticity of A and M, leads to an equality (Eq* — Eq) ({ijj\\M\\x} = 
for arbitrary two eigenstates and ||^)) belonging to eigenvalues E^ and 
Eq, respectively. This implies, in particular, that the energy eigenvalues Eq must 
be real as far as the eigenstate has non-zero norm ((x||A^||x)) ^ 0. Namely all 
solutions to our BS equation (2.25) for which the LHS of Eq.(2.28) is non-zero cor- 
respond to real eigenvalues Eq. In numerical work we actually find many solutions 
to (2.25) with complex eigenvalues Eq, but they all turn out to have zero-norm in 
accord with this observation. 

3. Decay Constant 

The pseudoscalar meson decay constant Fb is defined by 

%F B q^- iqx = (0| ^(x) lfll ^(x) |B(q)> \[W B , (3.1) 

where the factor of \[Mb is included so that the state |-B(q)) \[Mb satisfies the 
usual relativistic normalisation. Our definition would correspond to a pion decay 
constant F n = 93 x MeV. 

Setting r = x — y = in the definition of the BS amplitude (2.1) we find 

F B % = -N c [ tr[ 7At7 5x(p; q)] x • (3.2) 

Substituting the component expression of the BS amplitude x, we have 

F B ^fM B : = N c j *^L[A{u,x)-C{u,x)\, (3.3) 
which is an exact formula. 
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4. Isgur-Wise Function 



The expression for a heavy quark current operator ^'7^^ between two pseu- 
doscalar boundstates B' and B containing heavy quarks ^' and \1/ with masses M' 
and M may be described in terms of two form factors f(t) and g(t) as 



(5^)|*'(x) 7 ^0r)|5(q)> = [(„ + „%/(*) + (v - v'),g(t)} (4.1) 



where t — v ■ v'. When the initial and final mesons are the same (that is ^ = 
and Mb = Mb>), current conservation leads to g{t) = 0. The Isgur-Wise function 
£(t) W is defined as the leading term in the expansion of the first form factor f(t) 
in powers of m/M (or m/M') 



g(t) = (-)gi(t) + (-) g 2 (t) + -- - . 



The absence of the leading order term g$ (t) in g(t) can easily be seen : using 
the equality A + (tf' — ^)A + = 0, we find Xo(p'', — i>)xo(P'-, q) — (remember 

Xo = A+XO and xo = XoA+). 

£ (t) is independent of the heavy quark masses M and M' and is universal over 
pseudoscalar and vector heavy-light quarkonia: exactly the same function appears 
in the expansion of the form factors in the vector meson case. In the present 
approximation using constant mass for the heavy quark, the vector current matrix 
element (4.1) may be calculated from the diagram of figure 1. 
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1 

p'<9' 1 


— — ^ — ^ 


---\X(p'^)j 




q V J/ 






■p-r\q 



Fig.l Feynman diagram used to calculate the Isgur-Wise function. 



Thus we use 

(S'(q') I *'(0) 7/1 *(0) |S(q)> 

= Nc I^ tT ^ {P ' ] q ' MP] q)SLiP -^ q= P'- > (4 ' 3) 

where 77' = m/(M' + m) and £' = M'/(M' + m) for the final state meson £>'. 
Actually, in the equal mass case, the approximation of eq.(4.3) is consistent with 
current conservation to all orders: using i] + ( = 7/ + (' = 1, p — p' — rjq — rj'q' and 
M = M' we have 

q-q'=(v + ()q- W + C'W = (p + Cq) - (p' + CV) 

— ^-/ = (^ + a-M)-(/ + CV / -M)=5 H -^. (4 ' 4) 
Using this equation and the BS equation (2.13), it is easy to show that g(t) = 0: 
M B (v - v') 2 g(t) = (q- q'Yixll, ® S L || X » 

= ((x'\\(i-i')^s L \\x)) , 

(4.5) 

= «x'|| {Sr ® 5 L - 5^ ® S£) lx| (since S L = S{J 

= «x'll(^ - in«x» = o. 

We make one small comment: it is known that if a running heavy quark mass 
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function determined from the ladder SD equation were used, the matrix element 
(4.1) should be calculated by including gluon ladders also in the vector current 
channel, since otherwise current conservation would be violated. However dia- 
grams containing one or more gluon exchanges in this channel are easily seen to 
be suppressed by at least a power of {jfe) and so the approximation we use satis- 
fies current conservation to at least leading order in the m/M expansion. In fact 
go(t) = as we have seen before. 

The Isgur-Wise function is easily obtained by taking the leading order 
term of eq.(4.3) and contracting with (v + v')^\ we find 

2(1 + t)t(t) = N c f tr [Up'; q ) W + f)xo(p; q) s w (p -m = p- vW)] ■ (4.6) 
Jp' 

After taking the trace, we find the expression in terms of the invariant amplitudes: 
N c f du'x ,2 dx'd cos 9 f , - n ( L\\ ^12 \ / A \ 1 

where A' Q = Aq(u', x'), Aq = Aq(u, x) and so on, and the elements of the matrix L 
are given by 

L\\ — k\ + iu'(t + 1) 

L12 = -k\ - iu'ki(2t + 1) + (t + 1) [mki + iu'(t - l)(m - iv!) - x' 2 ] 

L21 = h(2m-iu) -x' 2 (t+ 1) (4.8) 

L22 = —k\{2m — iv!) + k\tx' 2 

2 2 

+ (£ + 1) [ki(iu' — m)(iu — 2m) + iu'x t — mx (t + 1)] 

with 



fci = • (y - f ; t) = pi v 7 ^ 2 - 1 • (4.9) 
We should note that the important condition £(i=l) = 1 is satisfied in our formu- 
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lation: actually, Eq.(4.6) may be rewritten into the form 



m = j^ix'oi^ ® s w \\ xo )) ■ 

where t becomes 1, v' and x'o become v and xo, respectively, so that the RHS of 
this clearly reduces to the normalisation condition (2.27) and gives 1. 

We evaluate the integral in (4.7) in the B' restframe in which v' = (1,0). In 
this frame we label the components of p' as 

/ /•// / /\ /•// /\ 

P = (l« ,Pl,P2,P3) = \ lu »Pl»P±) > 

(4.10) 

I'd I / I I'd 

Pi = x cose' , |p_|_| = x smb 1 , 

and take the integration contour of J d p'^ over the real axis of u' and p\. The 
momentum p of the initial state BS amplitude xip\ q) is determined by the relation 

p — rjq — p' — rj'q' — > p — p' — rj'q' + rjq . (4-H) 

and becomes complex. We require the corresponding arguments u and x for the 
component amplitudes Aq and Bq. Since u = p ■ v/i and x = a/— p 2 + (p ■ v ) 2 
are Lorentz invariants, it is simplest to calculate them in the w-restframe where 
v — (1,0). In this frame the vectors v' and p' take the form (since v • v' — t): 

v' = (t, v^T, 0, 0), 

(4.12) 

p' = (iu't + p'iVt 2 - 1, iu'y/t 2 - 1 + Pit, pj.) , 
so that p is given by 

p = p — r] q + r/q = p — mv +mv + (J{—) 

= (iut + p'i\/t 2 - 1 - m(t - 1), {iu -m)\/t 2 - l+p[t, pj_) + °(^) • 

(4.13) 

Therefore to leading order in m/M the arguments w and x of the initial state BS 
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amplitude Aq and Bo are given by 

= -p ■ v = ut + i[m(t — 1) — p'i \J t 2 — 1] , 



u 

i' 



x 



(4.14) 



2 = v/ [iuVi 2 - 1 + (pi* - mVt 2 - 1)] + p' ± 2 . 



Note that w and x in (4.14) are complex. In order to find the BS amplitudes 
Aq(u,x) and Bq(u,x) at these complex values, we need to perform an analytic 
continuation. Fortunately this continuation can be done using the BS equation 
itself: we simply put the complex values u and x into the BS equation (2.25) and 
find Aq(u, x) and Bq(u, x) by performing the RHS integration over v and y, which 
requires only the knowledge of the BS amplitude with real variables v and y. 



5. Heavy Quark Spin Symmetry of Leading Order BS Amplitudes 

The full BS equation eq.(2.13) of course applies to any heavy quark light anti- 
quark boundstate \B(q)). So the —1st order equation eq.(2.22) tells us that the 
leading order BS amplitude xo satisfies A_xo = for all boundstates (not just the 
pseudoscalar) ; xo oc A+. Therefore, for instance, the leading order BS amplitudes 
of pseudoscalar and vector boundstates generally have the following forms: 

XoPs = A+{A + B0) 15 , 

(5.1) 

Xov = A+ [{(Ay - Byf) + e • p{Cy + D v f)\ , 
where e p is the polarisation vector of the vector boundstate. 

Let us now show that our (improved) ladder BS equation satisfies the heavy 
quark spin symmetry in the leading order oim/M expansion, and that the pseu- 
doscalar and vector BS amplitudes are given by 

XoPs = A+l 5 {A-B]f) , 

(5.2) 

Xov = A+^(A - Bj) , 
with common functions A and B, and belong to a common energy eigenvalue. To 
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see this we rewrite the leading order BS equation eq.(2.24) in the form 

(p-v- E ) X o(pW - 2mA+) = / v^ X o(k)^g 2 C 2 D^(p - k) (5.3) 

Jk 

with xo subject to the constraint A_xo = 0. The important point is that in eq.(5.3) 
the heavy quark spinor index a of the amplitude xo a [3 (the left-leg) is left intact: 
there are no 7-matrix factors multiplying xo from the left. The leading order BS 
equation is therefore invariant under any heavy quark spin rotation which com- 
mutes with A_ — commutativity with A_ is required to maintain the constraint 
A-Xo = 0. In particular, in the restframe of the boundstate in which A_ = ^p^, 
it is invariant under 

xo(p;q) -> U(6)xo(p;q) , 

o-, (5.4) 
U{9) = exp (i0 ■ |) . 1 ; 

In this frame, the polarisation vector of the vector boundstate is of the form 
= (0,e) (|e| = 1) since e • v oc e • q = 0. Let us rotate the pseudoscalar amplitude 
through 6 = n around the eaxis by multiplying by U(ne). Noting the relations 

U (7rc) = ia e and A + o"75 = A + 7 , 

we find 

U{ne) A+75 (A - Bj) = iA+e ■ -y{A - Bp) . (5.5) 

But (aside from a constant overall factor —i) the RHS is nothing but the vector 
meson BS amplitude xov as given in eq.(5.2). So if the pseudoscalar BS amplitude 
XoPs i n eq.(5.2) satisfies the BS equation (5.3), then so does the vector amplitude 
Xov in (5-2), and the energy eigenvalues are the same. 

It is also important to note that the leading order BS equation (5.3) is totally 
independent of the heavy quark mass. This implies that the excitation energies 
determined as energy eigenvalues of that equation are also universal; namely, suffi- 
ciently heavy quark mesons with the same spin-parity or their partners under heavy 
quark spin symmetry should show almost the same excitation energy spectra. 
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In the same way as in the pseudoscalar and vector case, we can see that the 
scalar and axial vector mesons form a degenerate multiplet under the heavy quark 
spin symmetry, and that their leading order BS amplitudes are given in the form. 



Xos = A+(A S -B S ]£) , 



(5.6) 



6. Numerical Calculation 



6.1. Choice of Running Coupling Constant 

Asymptotic freedom requires that at scales /i 2 Aq CD the running coupling 
behaves (in the one loop approximation) as 



where the constant ao depends on the number of light quark flavours — we work 
with three light flavours for which ao — 47r/9. Unfortunately we have no guidance 
on the form of the coupling outside this deep Euclidean region. Previous investiga- 
tions of the improved ladder approximation 16 ' 101 have simply continued the form of 
(6.1) so as to have no singularities on the real axis in the spacelike region /i 2 > 0. 
However the prescriptions which were used are inadequate for the present calcula- 
tion of the Isgur-Wise function. This is because argument pj| + of ot{p\ + /cj|) 
becomes complex when we use the BS equation (2.25) to calculate x(u,x) at com- 
plex argument. We thus require that the coupling be analytic and that it have no 
singularities within the integration region used to perform the analytic continuation 
for the boosted BS amplitude. 



M^cd) 



(6.1) 
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We adopt the following form: 
a(/x ) - 



www ' ( , 2) 

Xq — X 



f(x) = x + k In 



1 + exp 



K 



with suitable real k > and xo > 1. A qcc j is a constant analogous to Aqcd which 
sets the scale of the coupling. Since we determine the value of A qc d mainly from 
the infra-red behaviour of the coupling, we should not expect it to agree with Aqcd 
which is usually fixed using data at much higher energies. Notice that f(x) ~ x 
for large positive x, so that the leading asymptotic behaviour of the coupling is 
correct regardless of the value of A qcc j. 

From eq.(6.2) one can show that a tends to the constant value a max = ao/ \ylxq 
in the low energy (or timelike) region x -C xo- We use the values 

x Q = 1.01, 1.05 and 1.10, 

corresponding to a m&x — 140.3, 28.6 and 14.6 respectively. 

Next we determine k. f(x) has branch point singularities at 

x s = xq + 27r/t(2n + 1) , n G Z , (6.3) 

with the branch cuts extending horizontally to the left if the principal value of 
the logarithm is taken; these singularities give rise to corresponding singularities 
in a(/i 2 ). So k must be sufficiently large that these singularities and branch cuts 
lie outside the integration region used to calculate the BS amplitude x( u , x ) with 
complex arguments u and x. From eq.(4.13), the complex momentum squared pj| 
is given by 



pi = u' 2 - 2im( t-l)u' + x' 2 - 2mx' cosO^ t 2 - 1 + 2m 2 ( t - 1 ) , (6.4) 
where u' and x' run over real values. We can avoid the singularity in / ( (p|, + 
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*D/ A qcd ) if > When Re (PE + fc E )/ A qcd = *0, 



|Im(p| + fc|)| = 2m(t - l)u' < 7r/tA 2 cd , for all u' ,x' > . (6.5) 

When Re(p| + A;|)/A2 cd = x , 

x A 2 qcd = Re (pi + fc|) 

= kl + u' 2 + (x' - m^t 2 - lcos^) 2 + m 2 ( 2t - 2 - (t 2 - 1) cos 2 9 ) 

> M '2_ m 2 (t _ 1)2) 



so that 



Thus we require 



(6.6) 

' < C = V / x A 2 cd + m2(t-l)2 . (6.7) 



2(t - 1) m / m 2 

for the range of values of t we consider. As will be seen below, at fixed k, m/A qcc j 
is determined by xq and the RHS is a function of only xq and t. Using the data 
for m/A qcc j obtained below for the three values of xq, we can estimate the RHS 
of (6.8) and find that k = 0.3 allows us to calculate the Isgur-Wise function for 
t < 1.55. As will be seen below, this region of t covers the kinematically allowed 
region for the semi-leptonic 5-decay which we are interested in. So we fix k equal 
to 0.3 henceforth. The reader may have noticed that there are also singularities in 
eq.(6.2) when f(x) = or f(x) = 1; it is easy to see that the condition (6.8) given 
above is sufficient to avoid these singularities. 
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6.2. Fixing the Mass Scale and the Light Quark Masses 



In order to treat dynamical chiral symmetry breaking consistently in the ladder 
approximation we should use the dynamical quark mass function obtained from 
the improved ladder approximation to the SD equation, and we should use the 
same kernel in both the SD and BS equations. It is easy to see that this approach 
correctly leads to massless pions in the chiral limit! 61 

In fact, the inclusion of the running mass is plagued with a number of technical 
difficulties related to the consistency of the axial WT identity 111 ' In the present 
calculation, there is the further problem that, after Wick rotation of the loop mo- 
mentum, the momenta flowing through the light and heavy fermion propagators, 
Snip + Co) an d Si,(p — r}q), become complex necessitating values of the quark mass 
function H(x) with complex arguments x = — (p + (q) 2 , — (p — VQ) 2 - These are 
not insurmountable problems (in particular the quark mass function with complex 
argument may be evaluated using the SD equation in the same way as we calcu- 
lated the boosted BS amplitude in Sect. 3), however they considerably complicate 
the formalism. In this paper we therefore use fixed quark masses and defer the 
discussion using the running mass function to a future paper. 

It turns out that the leading order calculations need no particular value of the 
heavy quark mass M, but we do have to fix the light quark mass m as well as 
the scale A qC( j used in the coupling. In order to do this we adopt the following 
procedure. First we solve numerically the improved ladder SD equation for the 
light quark mass function using the running coupling a{jp\ + k^) given in 



There is no unique definition of the constituent quark mass in terms of the quark 

[121 

mass function E; we work with the following two definitions 



eq.(6.2); 




(6.9) 




(6.10) 
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Next, to determine the energy scale A qcc j, we calculate the pion decay constant 
Ft^PS) using the obtained mass function E and the Pagels-Stokar formula 1131 which 
reads 

^( p s) = ^ y o dx — {x+nx?? — • (6-n) 

Since we know the Pagels-Stokar formula to agree with the ladder exact result to 
10 ~ 20%| 6,141 imposing the value i^(PS) = OSv^MeV allows us to fix A qcd . As 
a consistency check we also calculate the expectation value of the quark bilinear 
(ipip) 1GeV | 61 us i n g the formula 

9C 2 . 2 

/7,\ ( a ( A ) \ 11N *- 2N f N c f A , xE(x) 

with sufficiently large ultraviolet cutoff A, where Nf = 3 is the number of light 
flavours. This value (6.12) should be compared with the following "experimental" 
value given by Gasser and Leutwyler 



- (^)iGev = ( 225 ± 25 MeV ) 3 . 



The results obtained by this procedure are listed in Table. 1. 



x 


Aqcd 


mtype I 


mtype II 


(-<^> 1G ev) V3 


1.01 


638 


489 


288 


212 


1.05 


631 


482 


286 


213 


1.10 


625 


474 


285 


214 



Table. 1. Parameter values used in our calculation. Units are MeV. 
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6.3. Bethe-Salpeter equation 

In order to solve the BS equation (2.25) numerically, we discretise the variables 
u and x so that it reduces to a finite dimensional eigenvalue problem, which we solve 
using a standard linear algebra package. Since logarithmic scale seems natural, we 
discretise the variables U = ln(u/A qc( j) and X = ln(x/A qc d) at iVes = 30 points 
evenly spaced in the intervals 

U e [-10.0, 2.5] , X G [-4.5, 4.0] . (6.13) 

The integration kernel Ko(u,x;v,y), given explicitly in the appendix, has an in- 
tegrable logarithmic singularity at (u, x) = (v, y) which we avoide by using the 

[7] 

four point average prescription . Calculation of the binding energy Eq and the 
corresponding BS amplitude xo( u , x ) is straightforward, and we can immediately 
calculate the decay constant Fg by integrating the real part of Aq(u, x) via eq.(3.3). 
The number of sites Abs and the support of U and X in eq.(6.13) are large enough 
to guarantee the discretisation independence of the binding energy and the decay 
constant to within 1 % — see Table. 2. As a consistency check, we also confirm 
that for the fixed coupling hydrogen atom like case our program gives results which 
agree well with the formula Eq = ma 2 / 2 for small a. 



type 


XQ 


A^BS 


FbVMb 


p(0) 




I 


1.05 


18 


3468 


935.0 


498.3 


I 


1.05 


30 


3484 


931.7 


499.4 


I 


1.05 


34 


3486 


931.4 


499.3 



Table. 2. Abs independence. Eq and Eq are binding energies of ground state 
and first excited state, respectively. Units are (MeV) 3 / 2 for Fb\/Mb and MeV for 
4 0) and 

As can be seen in figure 2, the main supports of the integrands of normali- 
sation (2.28) and decay constant (3.3) are included in the ranges given in (6.13). 
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Unfortunately we can also see that these supports extend far into the infra-red - 
much further than in the case of the pion' 61 — and consequently our results depend 
on the infra-red behaviour of the coupling. We are therefore forced to include xq 
as a parameter in our model. 



Fig.2 



Our numerical results for a range of values of xq are shown in Table. 3. We 
prefer the value xq = 1.05 for the simple reason that at this value the shape of the 
running coupling most resembles the form used in our previous work*.' 6 ' 7,141 Ideally 
we would use the excitation energies — E^ to fix the parameter xq, but at 
present there is no experimental data for the masses of the excited pseudoscalar or 
vector B or D mesons. 



type 


XQ 


light quark mass 


EbVMb 


p(0) 






I 


1.01 


489 


2551 


1795 


1071 


724 


I 


1.05 


482 


3468 


935 


498 


437 


I 


1.10 


474 


4093 


666 


319 


347 


11 


1.01 


288 


2052 


1558 


972 


586 


11 


1.05 


286 


2738 


799 


447 


352 


11 


1.10 


285 


3205 


566 


279 


287 



Table. 3. xq dependence; Nbs = 18. Units are MeV except for Fb\[Mb which 
is in (MeV) 3 / 2 . 

The B and D meson decay constants and mass difference using a finer discreti- 
sation (iVes = 30) are given in Table. 4 where these meson masses are taken from 

* Indeed, we have calculated Fbs/Mb using the previous running coupling function 161 and 
confirmed that it produces almost the same result as the present one for the choice xq = 1.05. 
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Ref.[16]. 



type 


Xq 


mass 


^B(5279) 


^£>(1869) 


fi(lst) - 5(5279) 


I 


1.05 


482 


48.0 


80.6 


432 


n 


1.05 


286 


37.8 


63.4 


350 



Table. 4. Results for the meson decay constants and mass difference with Abs = 
30 and xq = 1.05. S(lst) denotes 1st excited state of pseudoscalar B meson. Units 
are MeV. 

As can be seen, our result for Fg is smaller than F %) and is much smaller than other 
values obtained using QCD sum rules, potential models and lattice simulations 

[17] 

(see the summary in the paper by Rosner ) — we note however that one gluon 

[18] 

exchange models tend to give small values, although even these are somewhat 
larger than our results. 

6.4. The Isgur-Wise function : £(t) 

After evaluating the BS amplitude of the ground state, we can now apply the 
formalism of section 4 and calculate the Isgur-Wise function £(t) by numerically 
evaluating the integral eq.(4.7). This integral itself is three dimensional, over u',x' 
and cos 6, but as explained in section 4, the calculation of the initial BS amplitudes 
A(u, x), B(u, x) with complex arguments u and x requires the evaluation of the BS 
equation (2.25) which involves a further two dimensional integral over v and y. So 
in fact we are calculating a five dimensional integral, which takes a lot of computer 
time when the number Abs of points on which the arguments u, x, v and y of the 
BS amplitudes are discretised becomes as large as 30. In order to save computer 
time, we evaluate the angle integral over cos# using the Gauss-Legendre formula 
for numerical integration which is known to give rather precise values with using 
a relatively small number of data points Aql- 

The main experimental interest at the moment is in semi-leptonic decay pro- 
cesses such as B -> D^lu, since these should allow calculation of the Kobayashi- 
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Maskawa matrix element V^b' 51 - For these processes, the invariant mass between / 
and v must be positive, and this bounds t from above 



t < 



M 2 B + M 2 r 



2M B M DM 



^-1.5-1.6. 



(6.14) 



We have chosen the parameter k in our running coupling so that we may evaluate 
the Isgur-Wise function within this region. The result can be seen in figure 3, 
where we have shown the Isgur-Wise function for the two cases of type I and type 
II masses. Figure 3 is the main numerical result in this paper. 



Fig.3 



In figures 4 and 5, we show how the result depends on the choice of the param- 
eter xq in the running coupling constant, for both the type I and type II masses. 
These data are based on the BS solutions on a coarser lattice with N^s = 18 points 
and cos 9 integration evaluated by Nq^ = 10 Gauss-Legendre formula. We see that 
the result is almost independent of xo for the type I case, whereas there is some 
dependence for the type II case. 



Fig.4 



Fig. 5 



Here a comment may be in order on some 'fluctuations' of our data points. We 
observe 'large' fluctuations occur for instance, at t — 1.5(xq = 1.01) point in Fig.4, 
t = 1.18(type I) point in Fig.3 and so on. We do not yet understand the precise 
origin of this phenomenon. But we suspect that some 'coherence' is occurring 
between the Gauss-Legendre formula and the discretisation for the BS amplitude 



2G 



data. Indeed, if we increase A<gl of the Gauss-Legendre formula, the Isgur-Wise 
function calculated by using the same BS data becomes more smooth as a whole, but 
the places at which such 'large' fluctuations occur change. Actually, for example, 
the large deviation at t — 1.5(xq = 1.01) in Fig.4 does not appear when Agl = 6. 
(Therefore the seeming 'oscillation' observed in Fig. 5 around t — 1.3 ~ 1.5 will 
also be a fake.) If we are allowed to average all the data obtained with various 
Agl(> 6), we obtain very smooth curves which coincide with the cited Agl = 10 
data with fluctuations smoothed out. 

The notable feature of these figures 4 and 5 is that the charge radius of the 
Isgur-Wise function 



is as large as 2.0 for type I mass or 1.8 for type II, and is quite insensitive to the 
choice of xq. These values of p 2 are certainly larger than the usual predictions of 
p 2 = 1.1 ~ 1.3. [17,19 ' 201 In our model the slope of the Isgur-Wise function does not 
decrease to 1.1 until region t — 1.1 ~ 1.3. If this feature of our result is correct, 
it would imply that V c b as extracted from the experimental data l21 ' 22 ' 23 ' 24 ' 2j ' 261 has 
been underestimated by several percent. 

To see this more explicitly, we fit our Isgur-Wise function £(£) of figure 3 to 

[27] 

the ARGUS data ( with setting their parameter tb = 1-32 ps ) by adjusting 
V ch so as to minimise X 2 = En ( \ v cb\ f (*n) ~ fnf where /„ and a n are the 
experimental data and standard deviations at t — t n . The result is shown in figure 



and see that the type I case gives a better overall fit to the ARGUS data. 



6. We find 




for type I case 
for type II case 



(6.16) 



Fig.6 
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Finally, we use our Isgur-Wise function with type I mass to estimate the p 

. f i r r 1 1 ,\ i r [21,22,24,28,29,30] 

parameters m various functional forms of £{t) proposed so tar. ine 
functional forms and the corresponding charge radii p 2 which give the best fit to 
our data are given in Table. 5, and the resulting curves are shown in figure 7. The 
model C gives the least x 2 , although the models B and C also give good fits. 



model 


function form 


P 2 


A 


l-p 2 (t-l) 


1.31 


B 


rfiexp [-(V-l)f^] 


2.16 


C 




2.03 


D 


exp[V(*-l)] 


1.87 



Table. 5. Various parametrisations of the Isgur-Wise function. The charge radii 
are extracted from our with N-qs = 30, xq = 1.05 and type I mass in Fig.3. 

Fig.7 
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APPENDIX : BS Kernel After Angle Integration 



After performing three-dimensional angle integration over cos0 = k ■ p/|k||p| 
the zeroth order kernel Kq(p, k) is found to be given by 

/ h-(u- vfh i(« - v)tfh - / 2 P ' k ) . 

where 

/- 1 1 1 ^(a; + y )2 + ( u _^2 

h = dcosO — — 79=7; — In 



(A.2) 



1 — (k — p) 2 2xy \ (x — y) 2 + (u — v ) 2 

1 1 2 



I2 = dcosO 



1 (A;-p) 4 ((x + ?/) 2 + (M-t;) 2 )((a;-y) 2 + (M-t;) 2 ) 

rP -k_ , fl k p x 2 + y 2 + (u - v) 2 

7f k = J^dcose^^= l -({x^ + ^ + {u-vf)h-h) . 



REFERENCES 

1. N. Isgur and M.B. Wise, Phys. Rev. Lett. 66 (1991), 1130. 

2. N. Isgur and M.B. Wise, Phys. Lett. B232 (1989), 113; Phys. Lett. B237 
(1990) 527. 

3. E. Eichten and B. Hill, Phys. Lett. 234 (1989), 511. 

4. H. Georgi, Phys. Lett. 240 (1990), 447. 

5. M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49 (1973), 652. 

6. K-I. Aoki, M. Bando, T. Kugo, M.G. Mitchard and H. Nakatani, Prog. 
Theor. Phys. 84 (1990), 683. 

7. K-I. Aoki, T. Kugo and M.G. Mitchard, Phys. Lett. B266 (1991), 467. 

29 



8. See for example N. Nakanishi, Prog. Theor. Phys. Suppl. 43 (1969), 1. 

9. K-I. Aoki, M. Bando, T. Kugo and M.G. Mitchard, Prog. Theor. Phys. 85 

(1991) , 355. 

10. P. Jain and H.J. Munczek, Phys. Rev. D44 (1991), 1873; Phys. Rev. D46 

(1992) 438. 

11. T. Kugo and M.G. Mitchard, Phys. Lett. B282 (1992), 162. 

12. H. Georgi and H.D. Politzer, Phys. Rev. D14 (1976), 1829. 

13. H. Pagels and S. Stokar, Phys. Rev. D20 (1979), 2947. 

14. T. Kugo and M.G. Mitchard, Phys. Lett. B286 (1992), 355. 

15. J. Gasser and H. Leutwyler, Phys. Rep. 87 (1982), 77. 

16. Particle Data Group, Phys. Rev. D45 (1992), SI. 

17. J.L. Rosner, Phys. Rev. D42 (1990), 3732. 

18. M. Suzuki, Phys. Lett. 162B (1985), 392. 

19. H.-Y. Jin, C.-S. Huang and Y.-B. Dai, Z. Phys. C56 (1992), 707. 

20. T. Mannel, W. Roberts and Z. Ryzak, Phys. Lett. B254 (1991), 274. 

21. P. Ball, preprint HD-THEP-92-25, to appear in Proceedings of "27th Ren- 
cotres de Moriond", 1992. 

22. J.G. K6rner and G.A. Schuler, Z. Phys. C38 (1988), 511. 

23. M. Bauer, B. Stech and M. Wibel, Z. Phys. C29 (1985), 637. 

24. B. Grinstein, N. Isgur, D. Scora and M.B. Wise, Phys. Rev. D39 (1989), 
799. 

25. K. Hagiwara, A.D. Martin and M.F. Wade, Nucl. Phys. B327 (1989), 569. 

26. J.M. Cline, G. Kramer and W.F. Palmer, Phys. Rev. D40 (1989), 793. 

27. ARGUS Collaboration, Z. Phys. C57 (1993), 533. 

30 



28. M. Neubert and V. Rieckert, Nucl. Phys. B382 (1992), 97. 

29. M. Neubert, Phys. Rev. D45 (1992), 2451. 

30. M. Neubert, Phys. Lett. B264 (1991), 455. 



FIGURE CAPTIONS 



Figure 2. Integrands of normalisation (2.28) and Fb (3.3). N-qs = 30 and 
xq = 1.05. The upper 9/10 of the figures are clipped. 

Figure 3. Isgur-Wise functions calculated with iVes = 30, Nql = 10 and 
xq = 1.05 for both type I and type E masses. 

Figure 4. xq dependence of the Isgur-Wise function with type I mass, N-q$ = 
18 and N GL = 10. 

Figure 5. xq dependence of the Isgur-Wise function with type E mass,iVBs = 
18 and N GL = 10. 

Figure 6. £(y) ■ |V^b| versus y. xq = 1.05 and both type I and II masses are 
used. 

Figure 7. Least squares fits of the four theoretical models A, B, C and D to 
our numerical data. 
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Fig.4 
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Fig.5 
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Fig.6 
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Fig.7 
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